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Preemptive  Scheduling  Of  A  Multiprocessor  System 

With  Memories  To  Minimize  L  * 

max 

Ten  Hwang  Lai  and  Sartaj  Sahni 
University  of  Minnesota 


Abstract 

We  develop  an  6( k  n  +  nlogn)  algorithm  to  obtain  a  preemp¬ 
tive  schedule  that  minimizes  Lmax  when  n  jobs  with  given 
memory  requirements  are  to  be  scheduled  on  m  processors 
Xn>m)'of  given  memory  sizes*  k  is  the  number  of  distinct  due 

dates.  The  value  of  the  minimum  L  can  itself  be  found  in 

max 

0(kn  +  nlogn)  time. 
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1 .  Introduction 


The  problem  of  scheduling  n  jobs  on  a  multiprocessor  system 
consisting  of  m  processors,  each  having  its  own  independent 
memory  of  size  J.  has  been  considered  by  Kafura  and  Shen 
[3].  Associated  with  each  job  is  a  processing  time  t^  and  a 
memory  requirement  itk  .  Job  j  can  be  processed  on  processor 
i  iff  nij  _<  J^.  No  job  can  be  simultaneously  processed  on 
two  different  processors  and  no  processor  can  process  more 
than  one  job  at  any  given  time  instance.  In  a  preemptive 
schedule,  it  is  possible  to  interrupt  the  processing  of  a 
job  and  resume  it  later  on  a  possibly  different  processor. 
In  a  nonpreemptive  schedule,  each  job  is  processed  without 
interruption  on  a  single  processor. 


Obtaining  minimum  finish  time  nonpreemptive  schedules 
is  NP-hard  even  when  m  =  2  and  =  J2  [2].  Hence,  Kafura 
and  Shen  [3]  study  the  effectiveness  of  several  heuristics 
for  nonpreemptive  scheduling.  For  the  preemptive  case,  they 
develop  an  O(nlogn)  algorithm  that  obtains  minimum  finish 
time  schedules  (without  loss  of  generality,  we  may  assume  n 
m) .  Their  algorithm  begins  by  first  computing  the  finish 
time,  f*,  of  a  minimum  finish  time  schedule.  This  is 

done  as  follows.  First,  the  jobs  and  processors  are  reor¬ 
dered  such  that  J,  >  and  m,  >  m0  >  •••  >  m  . 

This  reordering  takes  O(nlogn)  time  (again,  we  assume  n  >_ 
m)  .  Let  F^  be  the  set  of  all  jobs  that  can  be  processed 
only  on  processors  1,  2,  ...,  i  because  of  their  memory 

requirements.  Let  X^  be  the  sum  of  the  processing  require¬ 
ments  of  the  jobs  in  F, ,  X-  =  0  iff  F.  ^  Kufura  and  Shen 

[3]  show  that 

f*  =  max{  max^t^},  max^X-j  /  i } }  .  (1.1) 

The  jobs  may  now  be  scheduled  in  the  above  order 
(m^  2.  m2  —  1.  mn)  using  f*  and  McNaughton's  rule  [4]. 
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In  this  paper,  we  extend  the  work  of  [3]  to  the  case 
when  each  job  has  a  due  time  d^  associated  with  it.  Every 
job  is  released  at  a  common  release  time  c^.  We  are 
interested  in  first  determining  whether  or  not  the  n  jobs 
can  be  preemptively  scheduled  in  such  a  way  that  every  job 
completes  by  its  due  time.  A  schedule  that  has  this  pro¬ 
perty  is  called  a  feasible  schedule. 

In  Section  2,  we  show  that  the  existence  of  a  feasible 
schedule  can  be  determined  in  polynomial  time  using  network 
flow  techniques.  The  complexity  of  the  algorithm  that 
results  from  this  approach  is  0(kn( n+kr ) log  ( n+kr) )  where  k 
is  the  number  of  distinct  due  dates  and  r  the  number  of  dif¬ 
ferent  memory  sizes  in  {J^,  J2,  •♦*,  Jm) .  In  fact,  a  feasi¬ 

ble  schedule  (whenever  one  exists)  may  be  obtained  in  this 
much  time.  In  Section  3,  we  develop  another  algorithm  for 
this  problem.  This  algorithm  is  considerably  harder  to 
prove  correct  but  has  a  complexity  that  is  only  0(kn  + 
nlogn) .  A  feasible  schedule  can  be  constructed  in  0(k  n  + 
nlogn)  time.  In  arriving  at  the  algorithm  of  Section  3,  we 
develop  a  necessary  and  sufficient  condition  for  the 
existence  of  a  feasible  schedule.  With  the  help  of  this 
condition,  in  Section  4,  we  develop  an  algorithm  to  obtain  a 
schedule  that  minimizes  the  maximum  lateness.  This  algo- 
rithm  is  also  of  complexity  0(k  n  +  nlogn) . 

Sahni  [5]  and  Sahni  and  Cho  [6  and  7]  have  done  work 
related  to  that  reported  here.  They  have  considered  preemp¬ 
tive  scheduling  of  n  jobs  with  due  dates  when 

•  Jj  ■  • • •  *  Jm.  For  the  special  case  when  all  memory 
sizes  are  the  same,  Sahni  [5]  has  developed  an  O(nlogmn) 
algorithm  to  obtain  a  feasible  schedule  (when  one  exists) . 
Sahni  and  Cho  [6  and  7]  have  obtained  efficient  algorithms 
for  the  case  when  »  Jj  *  *  *  *  *  and  the  processors  run 
at  different  speeds. 
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2_.  The  Network  Flow  Algorithm 

In  this  section,  we  develop  a  conceptually  simple  algorithm 
to  determine  a  feasible  schedule  whenever  such  a  schedule 
exists*  This  algorithm  consists  of  two  phases.  Let 
c. ,  c0,  c,  be  the  distinct  due  times  in  the  multiset 
{d^  d2'  ***'  dn^ *  We  may  assume  that  c0  <  <  • * •  <  c^ 
(recall  that  c0  is  the  common  release  time).  In  the  first 
phase  of  our  algorithm,  we  determine  the  amount  6^  ^  of  job 
i  that  is  to  be  processed  in  the  interval  cj_i  to  cj*  This 
is  done  using  network  flows.  In  addition  to  the  source  (s) 
and  sink  (t)  nodes,  the  network  to  be  constructed  consists 
of  job  nodes,  interval  nodes,  and  summation  nodes. 

For  each  job  there  is  a  job  node.  The  source  node  con¬ 
nects  to  each  job  node  via  a  directed  edge.  The  capacity  of 
the  edge  to  job  node  i  is  t. .  If  job  i  has  a  due  time  Cj, 
then  there  are  exactly  j  interval  nodes  corresponding  to  it. 
These  nodes  represent  the  time  intervals 
[c0,cl],  [c^c 2]f  .  ..,  LCj^Cj].  From  each  job  node  there 
is  a  directed  edge  to  each  of  its  interval  nodes.  Each  such 
edge  has  a  capacity  equal  to  the  length  of  the  interval 
represented  by  the  interval  node.  In  Figure  2.1,  the  inter¬ 
val  node  labeled  (i,j)  corresponds  to  job  i  and  time  inter¬ 
val  [c , c j] . 

The  summation  nodes  are  divided  out  into  k  summation 
chains  with  each  chain  being  used  to  sum  up  the  flows 
through  each  of  the  k  sets  of  interval  nodes.  Let 
Q^,  ^2'  ***'  ^r  ^1  >  ^2  >  *  *  *  >  ^r^  the  distinct  memory 


sizes  in  the  multiset 

1^1  •  3  2  9 

•  •  •  ,  J  }  . 
mJ 

Let 

z^  be 

the 

number  of  processors 

of  size 

Q^,  and  let 

Pi  = 

i  Zj. 

3=1  3 

Let 

Qr+i  *  0.  There  are  exactly  r  summation  nodes  in  each  sum¬ 
mation  chain.  The  uth  node  in  chain  j  is  labeled  (u,j)  in 
Figure  2.1.  There  is  an  edge  with  capacity  Cj-Cj_^  from 
the  interval  node  (i,j)  to  the  summation  node  (u,j)  iff 


®u+l  <  mi  —  ^u'  From  no<^e  u  summation  chain  j  to  node 
u+1  there  is  an  edge  with  capacity  Pu(cj“cj_i)» 
l<u<r,  lili*.  The  last  node  of  each  chain  j  connects  to  the 
sink  node  and  has  capacity  m(  c  j-c  j_j_ )  . 


Job  Nodes  Interval  Nodes  Summation  Nodes 


Figure  2.1 


» rnmmmam 
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We  claim  that  there  is  a  feasible  schedule  for  the  n 

jobs  iff  the  maximum  flow  through  the  constructed  network  is 
n 

i  t.,  To  see  this,  first  observe  that  the  flow  cannot 
i31! 

n  n 

exceed  £  t. .  Suppose  that  there  is  a  flow  of  2  t.  in  the 
i=l  i=l 

network.  Let  6.  .be  the  flow  in  the  edge  from  job  node  i 
1  *  J 

to  interval  node  (i,j).  From  the  capacity  on  the  edges  from 

job  nodes  to  interval  nodes,  it  is  clear  that 

j  _<  cj-cj_i*  From  this,  the  capacities  on  the  edges  from 

summation  nodes,  and  Eq.(l.l),  it  is  clear  that  6.  .  ,  1  <  i 

1  •  3  ~ 

<_  n  can  be  scheduled  in  the  interval  [Cj_^,Cj]  using  the 

Kafura-Shen  algorithm,  1  £  j  <_  k.  Hence,  there  is  a  feasi- 

n 

ble  schedule  when  the  maximum  flow  equals  >  t..  it  is 

i=l 

readily  seen  that  the  reverse  is  also  true;  i.e.,  when  there 

n 

is  a  feasible  schedule  there  is  a  flow  of  value  £  t. . 

i=l  1 

So,  in  phase  1  of  our  algorithm  we  construct  the  flow 

network  and  determine  the  maximum  flow.  If  this  is  less 
n 

than  2  t . ,  then  no  feasible  schedule  exists.  If  the  max- 
i=l 

n 

imum  flow  equals  £  t. ,  then  in  phase  2  a  feasible  schedule 

i=l 

is  constructed  using  the  Kafura-Shen  algorithm  together  with 

the  interval  flows  6.  .  in  the  maximum  flow. 

1  *  J 


In  determining  the  complexity  of  this  approach,  we 
first  note  that  the  interval  nodes  are  easily  eliminated. 
Since  exactly  one  edge  enters  an  interval  node  and  exactly 
one  leaves,  these  two  edges  may  be  combined  into  one.  The 
total  number  of  nodes  in  the  resulting  network  is  n+kr+2 . 

The  number  of  edges  is  0(kn).  A  maximum  flow  in  a  v  node  e 
edge  network  can  be  found  in  O(evlog^v)  time  [1].  Hence, 
the  maximum  flow  in  our  network  can  be  determined  in 
0(kn(n+kr) log  (n+kr) )  time.  Following  this,  k  applications 
of  the  Kafura-Shen  algorithm  are  needed.  Hence  the  total 
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time  needed  to  find  a  feasible  schedule  (when  one  exists)  is 
0(kn(n+kr) log^(n+kr) ) . 

2 

3^  An  0(k  n  +  nlogn)  Algorithm 
2 

An  0(k  n  +  nlogn)  algorithm  (where  k  is  the  number  of  dis¬ 
tinct  due  times)  to  obtain  a  feasbile  preemptive  schedule 
for  a  set  of  n  jobs  may  be  arrived  at  in  the  following  way. 
As  before,  each  job  is  characterized  by  a  triple  (t^,d^,m^) 
where  t^  is  the  task  time  of  job  i;  d^  is  its  due  time;  and 
ra^  its  memory  requirement.  Let  c^,  C2»  •••«  c^, 

C1  K  c2  <  *  *  *  <  ck'  ^enote  the  distinct  due  times  in  the 
multiset  {d^,  d2»  •••»  dn) .  Let  Cq  be  the  common  release 
time  for  the  n  jobs.  Without  loss  of  generality,  we  may 
assume  that  c^  <  •  Further,  we  may  assure  that  the  jobs 
are  ordered  by  their  memory  requirements;  i.e., 
m.  >  m0  >  •  •  •  >  in  . 

The  m  processors  are  assumed  to  be  ordered  by  their 
memory  size,  i.e.,  l  J2  —  ' *  *  *  —  Jm*  Let 

Mg»0,  M^,  Mj,  •••»  Mr=m  be  such  that  processors 

M^+l,  M^+2,  have  the  same  memory  size  and 

JM  >  >  i  ,  0<i<r  (for  convenience,  assume  that  J_.1*0). 

The  processors  M^+l,  ...»  define  processor  class  i+l, 

0^i<r.  Similary,  let  ^=0,  ,  N2*  .  ..,  Nr=n  be  such  that 

jobs  N^+l,  N^+2,  ...»  have  a  memory  requirement  that  is 

larger  than  JM  , ,  but  not  larger  than  J..  ,  0<i<r.  It 

Mi+1+1  Mi+1  “ 

should  be  clear  that  a  job  j  such  that  <  j  _<  can  be 

processed  only  on  processors  1,  2,  ...,  .  Jobs 

N^+l,  ...»  define  job  class  i+l. 

It  is  easy  to  see  that  in  every  feasible  schedule  for 
the  n  jobs,  at  least 


b(i, j) 
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amount  of  job  i  must  be  completed  by  c^,  l£i£n,  0£j<k. 
Observe  that  if  there  exist  i  and  j  such  that  b( i, j ) >c j-Cg, 
then  there  is  no  feasible  schedule. 

Of  the  minimum  amount  b(i,j)  that  must  be  completed 
before  c^  ,  at  most 

a(  i,  j  )  =  min{b(i,  j)  .Cj^-Cg} 

can  be  completed  by  c^. 

Define  B(i,j)  to  be  the  sum  of  the  b(q,j)s  for  those 
jobs  q  in  job  class  i.  Define  A(i,j)  in  a  similar  manner. 
Specifically, 

Ni 

B(  i ,  j  )  *  1  b(q,  j  )  ,  l<i_<r,  0_<j<k, 

q=Ni_1+l 

and 

Hi 

A(i,  j)  =  I  a  ( q,  j  ) ,  l_<i_<r,  0£j_<k. 

q=Ni_1+1 

B(i,j)  gives  the  minimum  amount  of  job  class  i  that  must  be 
completed  by  c^.  A(i,j)  gives  the  maximum  amount  of  B(i,j) 
that  can  be  done  by  c^ . 

One  may  easily  verify  that 

0  b(i,j)  _<  b(i,j+l),  lj<i<jr»,  0_<j<k.  (3.1a) 

0  <_  a(i,j)  £  a(i,j+l),  l_<i<n,  0_<j<k.  (3.1b) 

0  ^  B(i,j)  <_  B(i, j+1),  l^i^r,  0^j<k.  (3.1c) 

0  <  A(i,j)  _<  A(i,j+1),  lj<i^r,  0^j<k. 


(3. Id) 


smm 


Lemma  X  obtains  an  interesting  inequality  between  A(i,j)  and 
B(i,j).  This  inequality  will  be  used  later. 


Lemma  1_:  If  t  <_  d^-Cg,  1  _<  q  £  n,  then 

<  liitli,  l <i<r ,  l<j<k. 


cj"cl 


'  Cj‘C0 


Proof;  Assume  that  t  _<  d^-Cg,  1  ±  q  _S  n.  Let  q  be  in 
the  range  [Ni_1+l,Ni].  Since  a(q,j)  =  min{b(q,j),  ci_c0}' 
a(q,j)  =  b(q,j)  or  a(q,j)  =  c^Cg. 

If  a(q, j )  ■  b(q,  j ) ,  then 


cl~c0 

cj-c0 


b  (  q ,  j  )  _<  a(q,  j  ) 


as  Cj. 


If  a(q,j)  =  ci_c0'  then  since  t  _<  d^-Cg 
b(q,j)  1  cj“c0'  we  9et 


implies 


c  —  c 

-b(q,  j)  <  c.-Cg  =  a(q,j) 

Cj-C0 


So,  in  both  cases  we  have 


Cl~C0t 

crc0 


•b(q,j)  <  a(q,  j  )  . 


Hence, 


:l‘c0 


Cj-c0q=N._1+l 


2  b(q, j )  <  2  a ( q , j ) 


q"“i-l+1 


cl_c0 


B( i, j)  <  A( i, j) 


or 
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ci“ c0  c-j  ca 

j)  1  ^L— ^A(i,  j) 

Cj  Cl  cj- ci 


3.~  0[B(i,  j)-A(i,  j)]  _<  B  (  i ,  j  ) 

Cj-Ci 


B  (  i  i  ]  )  A  (  i ,  ~] )  ^  B  (  i  ,  j  )  r  -j 
Cj”Cl  ”  Cj~C0 


Define  a  capacity  function  C(i,j)  such  that 


C(i,j)  -  )  ( c  j-c0  )  >  lj^i^r, 

C(  i*  j )  gives  the  available  processing  capacity  in  processor 
class  i  from  the  release  time  c0  to  the  due  time  c^. 

Let  n  be  the  set  of  all  nonincreasing  functions  O'  with 
domain  {0,  1,  2,  .  ..,  r}  and  range  {0,  1,  k}.  Recall 

that  r  is  the  number  of  processor  classes  and  k  the  number 
of  distinct  due  times.  Thus 


n  =  {  a  I  a:  {0, 1, . . . , r}-> {0, 1, . . . ,k}  and  o(  i)  >a(  i+1 ) ,  0_<i<r}  . 

ir  defines  the  set  of  profile  functions .  1  For  example,  con¬ 
sider  the  case  r  =  4,  k  =  5,  and  the  profile  function  <3  such 
that  or(0)  =  0T(1)  =  or(2)  =  4?  <3(3)  =  2;  and  or(4)  =  1.  Figure 
3.1  displays  <3  pictorially. 


One  readily  sees  that  if  the  given  job  set  has  a  feasi¬ 
ble  schedule,  then  it  must  be  the  case  that 
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Processor 

class  C0  C1  C2  c3  C4  C5 

1 

2 

3 

4 


Figure  ,3*i  Example  <3 

a  s 

I  B(i,cr(i))  <_  1  C(i,o,(i)) 

i=l  i=l 

for  all  s,  l£s£r  and  all  <3  <  ir.  In  particular,  if  there  is 
a  feasible  schedule,  then 

r  r 

1  B(i,<y(i))  <_  1  C(  i ,a(  i) )  (3.2) 

i=l  i=l 

for  every  <3  *  ir.  We  shall  show  in  Theorem  3.1  that  there  is 
a  feasible  schedule  iff  (3.2)  holds. 

We  note  that  for  any  a  *  it,  if  B(g,0”(q))  >  Clq.o'tq)) 

then  at  least  B(q,<y(q))  -  C(q,<y(q))  amount  of  class  q  jobs 
must  be  done  on  processor  classes  1,  2,  ...,  q-1  by  time 

cO(q)  ’  Hence  Btq^fq)  )  -  C(q,<y(q))  gives  a  lower  bound  on 
the  overflow  from  processor  class  q  to  processor  classes  1, 
2 ,  . . * ,  q- 1 . 
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Define  as  below: 

ir  j  =  {  a  I  o  <  n  and  o(i)_<j,  0^i£r}. 


We  see  that  for  any  profile  <3  * 


r  r 

2  B(q,o(q))  -  2  C(q,0'(q)  ) 

q=i+l  q=i+l 


gives  a  lower  bound  on  the  total  overflow  under  this  profile 
from  processor  classes  i+1,  .  ..,  r  to  processor  class  i, 
0<i<r.  Hence, 


r  r 

Z(i,j)  =  max  {  >  B(q,<r(q)  )  -  2  C(q,<r(q))}  (3.3) 

&*n.  q=i+l  q=i+l 

is  also  a  lower  bound  on  the  overflow  from  processor  classes 

r,  r-1  i+1  to  the  processor  class  i  up  to  time  c^. 

From  (3.3),  it  follows  that  Z(r,j)  =  0,  0^jjck.  Also, 
<  ^i-c0'  then  b(i,0)  =  0  ,  l_<i£n  and  we  obtain 


Z( i, 0)  =  0,  0<i<  r . 


(3.4) 


Furthermore,  since  <y:  {0, 1,  .  .  .  ,  r}->  {0}  is  in  for  all  j,  it 
follows  from  (3.3)  that  when  t^  <_  d^-c0,  l£i<n,  then 

Z(i,j)  _>  0,  0^ij^r,  0j<j_<k.  (3.5) 


From  (3.3)  we  may  obtain  a  sirtiple  recurrence  for 
Z(i,j).  Let  Of1  *  n j  be  the  3  at  which 

r  r 

5  B(q,<y(q))  -  £  C(q,cy(q)) 

q=i+l  q=i+l 

is  maximum.  Assume  that  i  <  r.  If  O'1  (i+1)  ^  j,  then  Z(i,j) 
*  Z(i,j-1).  If  O'*  (i+1)  =  j,  then 


-  13  - 


Z( i, j )-  2  B(q,a'(q))-  *  C(q,a' (q) ) 

q=i+l  q=i+l 

r  r 

=  2  B(q,0" (q) )—  2  C(q,a' (q) )+B(i+l, j)-C(i+l, j) 

q=i+2  q=i+2 


=Z ( i+1 , j ) +B ( i+1 , j ) -C ( i+1 , j  )  . 


This  yields: 


Z(i, j)= 


0 

Z(i,0) 

max{Z( i, j-1 ) ,  Z( i+1 , j )+B( i+1 , j ) -C( i+1 , j ) } 


if  i  *  r 
if  j  =  0 
otherwise. 


(3.6) 


Define  D(i,j)  as  below: 

J  0  j  =  0 

D(i»j)  -  |  c(i,l)  otherwise. 

Let  a  *  ir  be  a  profile  function.  B(q,<T(q)  )-A(q,<y(q)  ) 
is  a  lower  bound  on  the  amount  of  class  q  job  processing 
that  must  be  done  between  c^  and  ca(q)  •  Hence,  B(q,(y(q))- 
A(q,<y(q)  )-C(q,cr(q)  )+D(q,cy(q)  )  (abbreviated  as  [B-A- 
C+D]  (q,<y(q)  ) )  is  a  lower  bound  on  the  overflow  from  class  q 

jobs  from  time  c^  to  time  c^  ^  .  Consequently  X(i,j)  as 

defined  below: 


X( i, j)  =  max  {  2  [B-A-C+D](q,a(q) ) )  (3.7) 

a*  if  j  q=i+l 

is  a  lower  bound  on  the  overflow  from  processor  classes 
r, ...,i+l  to  the  class  i  from  time  c ^  to  Cj. 

Prom  (3.7),  it  follows  that 
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X(r,j)  =  0,  0£j<k. 


(3.8) 


Also,  if  _<  d^-c0,  l£i£n,  then  b(i,0)  =  a(i,0)  =  0,  and 

b(i,l)  =  a(i,l),  l<i<n.  From  this  and  (3.7),  we  obtain 


X( i, 0)  =  X( i, 1 )  =  0,  0<i<r . 


(3.9) 


From 


X(i, j) 


(3.7)  the  recurrence 

0 

=  X(i,0) 


i=r 

j=0  (3.10) 


|  max{X( i, j-1 ) , [X+B-A-C+D] ( i+1, j ) }  otherwise 

may  be  obtained  in  the  same  way  as  (3.6)  was  obtained  from 
(3.3). 


From  (3.10),  we  see  that  if  <  d^-c0,  l<i<n,  then 

X(  i,  j  )  _>  0,  0£i£r ,  0£j<k.  (3.11) 

The  difference  between  Z  and  X  will  play  an  important 
role  in  the  development  of  our  algorithm.  Define: 

Y(i,j)  =  Z(i,j)-X(i,j),  0£i£r,  0£j<k. 

We  establish  in  Lemma  3  an  important  inequality  on  Y. 
Before  proving  this  Lemma,  we  obtain  some  results  needed  in 
its  proof. 

Lemma  2  s  If  t  <  dq  -  c0,  1  £  q  £  n,  then 

(Cj-c0)X(i,  j)  <_  (Cj-Ci  )Z(  i,  j  ) ,  0£i£r,  0£j£k. 

Proof:  Assume  that  t_  <  d_  -  Cm  1  <  q  <  n.  The  proof  is 

-  q  —  q  0,  —  ^  r 

by  induction  on  i  and  j . 


A 


Induction  Based )  When  i  *  r,  X(i,j)  =  Z(if  j)  =  0. 

Induction  Hypothesis ( 1 )  Assume  that  the  inequality  is 
correct  for  all  j  when  0  <  i  =  p  r. 

Induction  Step( 1 )  When  i  =  p-1,  the  inequality  may  be 
shown  correct  by  induction  on  j . 

I^.B.  {2)  When  j  =  0  or  1,  X(p-lf  j)  =  0  and  Z(p-l,j)  >_  0* 

I^.H.  (2^)  Assume  that  the  inequality  is  correct  when  k>j=q>l. 

1.8.  (2)  Let  j  =»  q+1 .  From  Eq.  (3.10),  we  see  that  there 
are  two  possibilities  for  X(p-l,q+l). 

Cased)  X(p-l,q+l)  =  X(p-l,q):  In  this  case, 

(Cq“c0)X(p-l,q+l )  =  (cq-c0)X(p-l,q) 

-  (cq"ci )Z(p-l,q)  (I .H. (2 ) ) 

1  (cq-Cl)Z(p-l,q+l)  (Eq .(3.6)) 

Since  Cq  -  c0  >  cq  -  c^  0,  we  get 

X (p-1 ,  q+1 )  _<  Z  ( p-1 ,  q+1 ) 

or 

*cq+l_cq)X*p_1'q+1 *  -  <cq+l“cq)z(P“l*q+l) • 

Adding  this  inequality  to  the  previous  one  yeilds: 


<  Z(p-l,q+l  ) 
Cq+1  ~  C0 


(Eq . (3.10))  [] 


Corollary  1:  If  t  <_  dq-c0'  1  1  *3  1  n»  then 

Y(i,j)  ^0/  0^i£r,  0£j£k. 

Proof:  For  j  >_  1,  this  follows  from  Lemma  2  and  the  fact 
Cj-c0  >  Cj-c^  _>  0.  For  j  =  0,  this  follows  from  Eqs .  (3.4) 
and  (3.9).  [ ] 


Lemma  3:  If  t_  <  d_  -  c„  1  <  q  <  n,  then 
-  —  q  —  q  0,  —  ’  - 

Y(i,j)  >_  Y(  i ,  j-1 ) ,  0^i£r ,  l^jfk. 

Proof:  Assume  that  t  _<  dq  -  c0,  1  _<  q  <  n.  The  proof  is 

by  induction  on  i  and  j . 

I  .B.  (_1 )  When  i  «  r,  Y(i,j)  =  Y(i, j-1)  =  0,  l^j<k. 

I  .H.  (_1 )  Assume  that  Y(i,j)  j>  Y(i,j-1)  for  all  j,  when  1  <_  u 
<  i  <  r  where  u  is  an  arbitrary  integer  in  the  range  [l,r- 
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13. 

1.5. (l)  We  need  to  show  that  Y(u,j)  Y(u,j-1),  l£j<k. 
I^B.(£)  When  j=l,  Y(u,l)>0  (Corollary  1)  and  Y(u,0)=0. 

I^.H.  (2^)  Assume  that  Y(u,j)  Y(u,j-1),  l_<j<b. 

1.5. (2)  We  need  to  show  that  Y(u,b)  >_  Y(u,b-1). 

Case  (i)  Z(u,b)  _>  Z(u,b-1)  and  X(u,b)  =  X(u,b-1): 

In  this  case,  it  is  readily  seen  that  Y(u,b)  _>  Y(u,b-1). 

Case  (ii)  Z(u,b)  =  Z(u,b-1)  and  X(u,b-1)  <  X(u,b)  » 

[X+B-A-C+D ] ( u+1 , b ) : 

This  case  is  not  possible.  To  see  this,  suppose  that  this 
case  is  possible.  Let  a  u  be  the  largest  a  for  which 

Z(a, j)=Z(a, j-1)  and  X(a, j-1 ) <[X+B-A-C+D] (a+1, j )  (3.12) 

for  some  j.  Let  c  be  the  smallest  j  for  which  (3.12) 
holds.  Note  that  c  >  0.  So, 

Z(a,c)aZ(a,c-l )  and  X( a, c-1 ) <[X+B-A-C+D] ( a+1 , c) .  (3.13) 

Since  X(a,c-1)  _>  0,  it  follows  from  Eqs .  (3.10)  and 

(3.13)  that  X(a,c)  >  0.  Assume  that 

Z(a,c)  =  Z(a,c-1 )  =  Z(a,c-2)  =  •••  =  Z(a,v)  ^  Z(a,v-1). 

Then,  it  follows  from  our  choice  of  a  and  c  that 

X(a, c-1 )  =  X( a, c-2 )  «  •••  =  X(a,v). 


If  v 


0,  then  0  =*  Z(a,v)  =  Z(a,c)  _>  X(a,c)  >  0.  Hence,  v  ^ 


jjgji  %  i  I  f 


0.  Now,  from  the  choice  of  v,  we  get 


Y(a,v)  =  Z(a,v)  -  X(a,v) 


<  [Z+B-C] (a+1, v)  -  [X+B-A-C+D] (a+1, v) 


=  [Y+A-D] (a+l,v) 


(3.14) 


Also, 


X( a, c-1 )  <  [X+B-C-A+D] ( a+1 , c ) , 


Z(a,c-1 )  >  [Z+B-C] (a+1, c) 


Y( a, v)  =  Y ( a, c-1 )  >  [Y+A-D] ( a+1 , c ) 


Substituting  into  (3.14)  yields: 


[ Y+A-D ]( a+1, v)  >  [ Y+A-D ]( a+1, c) 


Y(a+1, v)  >  [ Y+A-D ]( a+1, c)  -  [A-D](a+l,v) 


>.  Y(a+l,c)  (Eq.  (3.1)  and  definition  of  D) 


But,  a  +  1  >  u  and  so  from  I.H.(l),  it  follows  that 


Y(a+l,c)  >  Y ( a+1 , c-1 )  >  •••  >  Y(a+l,v). 


So,  case  (ii)  is  not  possible. 
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Case  ( iii)  Z(u,b)  =  [Z+B-C] ( u+1 , b)  and  X(u,b)  »  [X+B- 

A-C+D] (u+1 ,b) : 

Now,  Y(u,b)  =  [Y+A-D3 ( u+1 , b) .  Suppose  that 

Z(u,b-1 )  =  Z(u,b-2)  =  ...  =  Z(u,v)  #  Z(u, v-1 ) .  (3.15) 

From  the  proof  of  case  (ii),  it  follows  that  X(u,b-1)  * 

X( u,b-2 )  *  ...  «  X(u,v) .  So,  Y(u, b-1 )  =  Y(u,v).  If  v  =  0, 

then  Y(u,b-1)  =  Y(u,0)  =  0  Y(u,b).  If  v  ^  0,  then 

Y(u,b-1 )  =  Y(u,v) 

<_  CY+A-D]  (u+1  ,  v)  (Eqs.  3.15,  3.6,  3.10) 

_<  Y(u+l,b)+[A-D](u+l,v)  (l.H.(l)) 

±  [Y+A-D]  (u+1  ,b)  (Eq.  3.1) 

=  Y(u,b) .  [ ] 


We  are  now  ready  to  describe  our  preemptive  scheduling 
algorithm.  The  jobs  will  be  scheduled  in  k  phases.  In 
phase  j  we  determine  the  amount  of  each  job  i  that  is  to  be 
scheduled  from  to  Cj.  Once  this  amount  has  been  deter¬ 
mined,  the  actual  schedule  from  cj_i  to  Cj  is  constructed 
using  the  Kafura-Shen  algorithm. 

Procedure  COMPUTE_W  determines  the  amount  w^  of  job  i 
that  is  to  be  scheduled  from  Cg  to  c^ .  In  this  procedure, 
R0  denotes  the  amount  of  idle  time  remaining  on  processor 
classes  1,  2,  ...,  s  following  the  scheduling  of  the  w^s 
corresponding  to  jobs  in  job  classes  1,  2,  ...,  s.  Lemma  4 
obtains  some  interesting  properites  of  R  .  When  actually 
implementing  COMPUTE_W,  the  subscripts  on  h,  R,  and  Q  may  be 
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line 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


11 

12 

13 

14 

15 


Procedure  COMPUTE_W 

//w^  is  the  amount  of  job  i  to  be  processed  from 

c0  to  cl// 

R0  <"  0 

for  s  <-  1  to  r  do  //consider  jobs  by  classes// 
Qs  <-  (  j  I  A(s,j)  +  Y(s,j)  £  Rs-1  +  C ( s , 1 ) } 
if  Qg  *  (6  then  print(  '  infeasible  job  set') 
stop  endif 

h  <-  max{  j  I  j  *  Q  } 


s 

case 


:(D  hg  =  k;  <-  a(i,k),  Ns-1  <  i  _<  Ng 
:(2)  h  <k,  A(s,h  +1)+Y(s,h  )>R  ,+C(s,l): 

set  w^  such  that 

a(  i,hs)— wi— a(  i'hs+l )  •  Ns-l<;LlNs  and 


N. 


i  wi+Y(s,hs)  =  Rs-i+C(s,l) 


i(3)  else:  Wj  <-  a(i,  hs  +  1),  Hg.].<iiNs 
end  case 

N 


"s  "  *s-l  +  C<s'1>  -  „  £  .."i 

end  for 
end  COMPUTE  W 


Ns-1+1 


Figure  3.2 


omitted.  We  have  kept  them  in  the  version  given  in  Figure 
3.2  so  that  we  may  easily  refer  to  the  values  of  h,  R,  and  Q 
during  different  iterations  of  the  for  loop.  One  should 
also  note  that  in  case  (2),  since  A(s,h  +1)  +  Y(s,h_)  > 
R8_l  +  C(s,l)  and  A(s,hg)  +  Y(s,hs)  £  Rs-1  +  there 
exist  w.  ,  a(i,h_)  <  w.  <  a(i,h  +1),  such  that 

X  S  "  X  *“  s 
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These  w^s  are  easily  determined  by  first  setting  all  w^  *» 
a(i,hg),  Ng_^<i£Ng  and  then  incrementing  the  w^s  one  by  one 
(up  to  at  most  a(i,hg+l))  until  the  desired  equality  is 
satisfied . 


r  r 

Lemma  4:  If  >  Bti/O^i))^  j>  C (  i , <r(  i )  )  for  every  O'** it,  then 

i=l  i=l 

(1)  Qg  4  6  and  hg  £  1,  1  £  s  £  r. 

(2)  Rs  hs^'  1  1  s  1  r* 

(3)  R  <  Y(s,h  +1)  if  h  4  k,  1  <  s  <  r. 

S  *'■  S  5  “  ■ 


Proof:  Assume  that 


r  r 

I  B( i ,cr( i) )  £  I  c(i,<r(i) )  (3.16) 

i=l  i=l 


for  every  Cf  *  it.  In  particular,  using  a(i)  =  0,  1  £  i  £  r 

r 

in  (3.16),  we  obtain  I  B(i,0)  £  0.  Since  B(i,0)  £  0,  l£i 

i=l 

£  r,  it  follows  that  B(i,0)  =  0,  1  £  i  £  r  and  so  b(i,0)  =» 

0,  t^  £  d^-Cg,  and  b(i,l)  £  c±-cq,  1  £  i  £  n.  Hence, 

A(i,l)  =  B(i,l),  1  £  i  £  r.  Prom  this  and  Eq.  (3.9)  we 

obtain 


A( i , 1 )+Y ( i , 1 )  =  B(i,l)+Z(i,l),  l<i<r.  (3.17) 


Now  we  shall  show  (1),  (2),  and  (3)  by  induction  on  s. 

^.B.  From  (3.3)  and  (3.16),  it  follows  that  Z(0,1)  »  0. 
From  (3.6),  it  follows  that  Z(l,l)  +  B(l,l)  -  C(l,l)  £ 
Z(0,1)  £  0.  Combining  this  with  (3.17),  we  obtain  A(l,l)  + 
Y(l,l)  £  C(l,l).  Since  R0  =  0,  it  follows  that  1  4  Q  . 

Hence,  Q ±  t  6  and  hj  >  1* 
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.  i 


Ni 


By  definition,  Ri  *  R0  +  °( 1, 1 )  -  >  From  the 


V1 

N1 

algorithm  we  see  that  if  h,  =  k,  then  I  w.  =  A(l,h-.). 

V1 

But,  from  line  3,  it  follows  that  Y(l,h^)  R0  +  C(l,l) 

A(l,h.) .  Hence,  R^Yfl.hj).  If  hx<k  and  AU.^+l)  +  Y(l, 
h.)>_  Rp.  +  C(l,l),  then  from  line  10,  we  obtain  Y(l,h,)  = 

nl 

R0  +  C(l.l) 


I  wi.  Hence,  R^  Y(l,hj_).  Since 


N0+1 


wi  <_  a(i,h^+l)  in  this  case,  2  w^<A( i ,h^+l ) .  Also,  from  the 
definition  of  h^,  it  follows  that  A(  1,  hj_+l  )+Y(l  ,h]_+l )  > 


R0+C (1,1).  Hence,  R±  <  Y(l,h1+1) 


In  the  third  case  of  the  algorithm,  A( 1 ,h^+l )+Y( 1 ,h^ )  < 
R0+C (1,1)  and  Iwi  =  All.h^l).  Hence,  Yfl,^)  <  R0  +  C(l,l) 
-  S  w^  and  R^>Y(l,h^).  From  the  definition  of  h^  and  w^,  it 
follows  that  R^<Y(l,h^+l)  (same  proof  as  for  previous  case). 


I..H.  Assume  that  (1),  (2),  and  (3)  are  true  for  some  arbi¬ 
trary  s,  s  =  q,  1  _<  q  <  r. 


jl.S.  We  shall  show  that  (1),  (2),  and  (3)  are  true  when  s  = 
q+1 . 

R^  Y(q,hq)  (induction  hypothesis) 

>_  Y(q,  1 )  (Lemma  3  ) 

*  Z(q,l)  (X(q,l)=0  from  Eqs .  3.16  and  3.9) 

_>  [Z+B-C]  (q+1, 1 )  (Eq.  3.6) 

»  [Y+A-C] (q+1 , 1 ) .  (Eq.  (3.17)) 

Hence,  1  •«  Qq+i  and  so  Qq+^  ^  £  and  hq+^  >_  1 .  (2)  and  (3) 
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can  now  be  proved  in  the  same  way  as  in  the  induction  base. 

[] 

Before  establishing  the  correctness  of  our  scheme  to 
compute  the  w^s,  we  obtain  some  relationships  concerning  the 
amount  of  processing  t' ,  of  job  i  that  remains  to  be  done 
following  time  c.  .  Note  that  t*.  =  t4  -  w4  1  <  i  <  n. 

Define  b'(i,j),  a'(i,j),  B'(i,j),  and  A'(i,j)  to  be  the 

values  obtained  for  b,  a,  B,  and  A  when  tV  is  used  in  place 

of  t..  Let  C'(i,j)  =  (M.-Mj  i  )(cj-c,  ),  l_<i<^r,  l<j<k,  and 

Nt  J 

W(i)  *  I  w  ,  l<i<r.  Lemmas  5  and  6  enable  us  to  estab- 
q-Ni_1+l  <* 

lish  Theorem  1 . 

Lemma  5^  If  the  condition  of  lemma  4  is  satisfied  then 

(1)  B 1  ( i,  j )  <_  B( i, j )-A( i, j ) ,  j  _<  hi. 

(2)  B‘ (i,  j)  *  B(i,j)-W(i),  j  >  hi. 

Proof:  It  is  easy  to  see  that  for  any  job  q,  1  _£  q  ±  n, 

b’(q,j)  =  max{0,  b(q,j)-wq}. 

When  j  <_  h^,  a(q,j)  £  a(q,h^)  (from  Eq.  (3.1)).  Prom  the 

algorithm,  we  see  that  w^afq,!^)  in  cases  (1)  and  (2)  and 

wq  m  a(q,h^+l)  in  case  (3).  Hence,  in  all  cases,  we  have 
wq  >a  ( ,  j  ) ,  Ni-1+1  1  q  <  Ni,  j  <  For  B'(i,j),  j  <  hi( 

we  now  obtain 

B* (i, j)  -lb1 (q, j) 


imax(0,  b{q,j)-w 1 
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£  £max{0,  b(q,j)-a(q,j)} 

=  E(b(q, j ) -a(q, j ) )  (as  b(q, j )>a(q, j )^0 ) 

=  B( i, j  )  -  A( i, j ) . 


We  now  consider  the  case  j  >  ^  k.  From  cases  (2) 
and  (3)  of  the  algorithm,  Eq  (3.1),  and  the  definition  of 
a(q,j),  we  see  that 


a(q,h^)  £  Wq  £  a(q,h^+l)  £  a(q,j)  _<  b(q,j). 


So, 


b'(q,j)  *  max (0,  b(q,j)-w^}  =  b(q,j)-wq. 


Hence, 


B' (i, j)  =  B(i, j)-W(i) .  C  3 


For  convenience  in  proving  the  next  lemma,  we  define 
h0  =  1. 

r  r 

Lemma  6:  If  2  B(  i ,  0(  i)  )  <_  2  C(i.O(i))  for  every  O  <  it,  then 
i=l  i-1 

the  following  are  true  for  every  s,  0  ±  s  £  r  and  every  O 
such  that  O’  ( s )  >  1  : 


s  s 

(1)  I  B 1  ( i,o(  i) )  +  Z(  s,o(  s)  )  -  Rg  <_  2  C '  (  i/0(  i)  ) , 

i=l  1=1 


s  s 

(2)  If  0(  s)  <h  ,  then  1  B '  (  i,o(  i)  )+X  (  s.O(  s)  )  j<  I  C’(i,o(i)). 
"  3  i=l  i=l 


Proofs  Assume  that 
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r  r 

I  B(  i ,  0”(  i)  )j<  I  C(i,<3’(i))  for  every  <7  *  w.  (3.18) 
i=l  i=l 


We  shall  use  induction  on  s. 


I.B.  Prom  (3.3),  we  obtain 


r  r 

Z(0,a(0) )  =  max  II  B(  i,a‘  (  i)  )-  l  C(  i  ,a'  ( i)  )  K 3 . 19 ) 

°"4"c.(0)  i=1  1=1 


Using  <7'  such  that  o' (i)=0,  0  <  i  <  r  in  (3.18)  yields 
B(i,0)  =  0  and  hence  t  _<  d^-Cg,  1  <_  q  ±  n.  From  (3.5), 
(3.18)  and  (3.19),  we  now  get  Z(0,O'(0))  =  0.  From  Corollary 
1  and  (3.11),  we  get  0  =  Z(0,<T(0))  X(0,<y(0))  >_  0.  Hence, 

X(0,  <y(0))  =  0. 


0  0 

Clearly,  I  B'(i,cr(i))  =  I  C’(i,o,(i))  =  R_  =  0. 

i=l  i=l  v 


Hence,  when  s  =  0, 


s  s 

I  B' (i,<y(i)  )+z(s,<y(s)  )-R  <  I  C'(i,cr(i)) 

i=l  s  “  i=l 


and 


s  s 

I  B‘  (i,<y(i)  )+X(s,<y(s) )  <_  I  c '  ( i,<y(  i) ) 
i=l  i=l 


for  every  <7  ■*  n . 

I^.H.  Assume  that  (1)  and  (2)  are  true  for  s  =  t,  where  t  is 
in  the  range  0  _<  t  <  r. 

I^.S.  We  proceed  to  establish  (1)  and  (2)  when  s  =  t+1  by 

considering  the  three  cases  cr(t+l)  >  h  cy(t+l)  <_  h.+,  and 

<7( t)  hfc  ;  and  <y(t+l)  _<  hfc+1  and  <?(t)  >  hfc. 
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Case  Is  <y(t+l)  >  hfc+^:  We  first  obtain  the  following 
t  t 

I  B'U.crU))  +  Z(t,o-(t))  -  R.  <  1  C*  (i,<y(i) )  (i.H.) 

i=l  ”  i=l 

Rt  “  Rt+1  =  w(t+1)  ~  C(t+l,l)  (def.  of  Rt+1) 

B  '  ( t+1 ,  <y(  t+1 )  )  +  W(  t+l )  =  B  ( t+1 ,  <y(  t+1 )  )  (Lenuna  5) 

z( t+i , <y( t+l )  )+B(t+i,<y( t+l) ) -c ( t+l , <y( t+l ) ) 

£Z(  t,cr(t+l )  )  (Eq  .(3.6)  and  CT(t+l)  >1 ) 

<z(t,<y(t))  (Of(t)  >a(t+l )  and  Eq .(3.6)) 


Adding  these  four  equalities  and  inequalities  yields; 


t+1  t+1 

2  B'(i,or(i))%+  z(t+l,cr(  t+l ) )  -  R  ,  <_  I  C*  (i,<y(i) ) 
i=l  i=l 


Case  2:  <y(t+l)  _<  hfc+^  and  cJ(t)  <_  hfc ;  From  the  induction 

hypothesis,  we  have 

t  t 

I  B'(i,<y(i))+X(t,or(t))  <  2  C1  (i,cr(i)  )  .  (3.20)' 

i=l  i=l 

From  Eq.(3.10)  and  the  fact  that  cy(t)  a(t+l)  ^  1,  we  get: 

x(t,<y(t))  _>  x(t,cr(t+l ) )  >_  [X+B-A-C+D](t+I,a(t+1)). 

Using  Lemma  5,  this  reduces  to 

x(t,<y(t))  _>  x(t+i,cy(t+i)  )+b‘ ( t+i,<y(t+i )  )-c 1  (t+i,<y(t+i) ) . 

Combining  with  (3.20)  yields: 
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t+1  t+1 

I  B‘ (i,<y(i)  )+X(t+l,<r(t+l)  )  <  I  c  '  ( i ,  <7(  i )  )  .  (3.21) 

i=l  "  i=l 

Since  >_  Y(t+l,ht+^)  >_  Y( t+1 , <y( t+1 ) )  (Lemmas  4  and  3) 

we  conclude  that  Z(t+1,  <y(t+l))  -  Rfc+1  <  x( t+l,<y( t+1 ) )  . 

Substituting  into  (3.21)  yields 

t+1  t+1 

l  B’ (i,a(i) )+Z(t+l,a(t+l) )-R.  <  5  C* (i,a(  i)  )  . 

i=l  t+1  “  i=l 

case  1*  <J(t+l)  <_  hfc+1  and  <j(t)  >  hfc:  From  the  induction 

hypothesis,  we  get 

t  t 

1  B*  (i,a(i)  )+Z(t,<j(t)  )  -  R.  <  2  C'  (i,cr(i)  )  .  (3.22) 

i=l  “  i=l 

Since  hfc  ^  k,  we  obtain  from  Lemma  4: 

Rt  i  Y(t,ht+1)  <  Y(t,<y(t)  )  .  (3.23) 

From  Lemma  5,  Eq.  (3.10),  and  the  inequality  &(t)  >_  o(t+l)  > 
1 ,  we  get 

[X+B’-C'  ](t+l,a(t+l) )  £  X(t,cr(t+1))  _<  X(t,<y(t)  ) .  (3.24) 

Adding  (3.22),  (3.23),  and  (3.24)  yields: 

t+1  t+1 

2  B '  ( i ,  <y(  i )  )+X  ( t+1 ,  a{  t+1 ) )  _<  2  C'(i,<r(i)). 

i“l  i=l 

Using  the  same  reasoning  as  in  case  2,  we  may  now  conclude 
the  truth  of  (1)  when  s  =  t+1.  [] 

Theorem  1^  There  exists  a  feasible  preemptive  schedule  for 
the  given  n  jobs  if  and  only  if 

r  r 

I  B(i,<y(i))  <_  2  C(i,<y(i))  for  every  a  +  w  . 

i-1  i-1  K 
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Proof:  We  have  already  pointed  out  that  if  a  feasible 

schedule  exists,  then  the  above  inequality  is  satisfied  for 
every  cT  <  tt^.  So,  we  need  only  show  that  when  the  above 
inequality  is  satisfied  for  every  O'  *  n^,  there  is  a  feasi¬ 
ble  schedule.  Assume  that 

r  r 

2  B( i ,  o(  i) )  <  2  C( i,a(  i) )  for  every  a  *  (3.25) 

i=l  i=l 


Prom  (3.25)  it  is  clear  that  when  k  -  1,  the  t^s  and  Eq 
(1.1)  yield  f*  _<  c^-c^  and  so  a  feasible  schedule  exists. 


For  the  induction  hypothesis,  we  assume  that  there 

exists  a  feasible  schedule  when  (3.25)  is  satisfied  and  k=*q 

for  some  q,  l<c[.  We  show  that  if  (3.25)  is  satisfied  when 

k=*q+l,  then  there  is  a  feasible  schedule.  From  Lemma  4,  we 

see  that  Q  ^  ^  for  any  s .  Hence  procedure  COMPUTE_W  suc- 
s 

cessfully  computes  the  w.  s.  It  is  clear  from  COMPUTE_W  that 

Ns 


w..  <^a(  i,u)^c^-C0  where  u  =  hc+l 


i— 


or 


q+1  and  that 


2  W.  < 

i=l  “ 


M  (c,-ca)  -  R  <  M  (c,-c_)  for  every  s.  Hence,  the  w • s 

S  1  v  3  —  SlfcJ  -L 

satisfy  (1.1)  (i.e.,  f*  £  c^  -  Cg )  and  may  be  scheduled  from 
Cg  to  c^  using  the  Kafura-Shen  algorithm. 


Now,  consider  the  t1 .s.  We  know  that  X(r,j)  =  Z(r,j) 
■  0,  0  <  j  <  q+1.  If  hr<q+l,  then  from  Lemma  4,  we  obtain  0 
<_  Rr  £  Y(r,hr+1)  =  0  or  Rr  =  0.  Using  this  in  Lemma  6 
yields : 


r  r 

2  B'  (i,<y(i) )  <  S  C'  (i,df(  i)  )  for  every  <5  ■«  w_+1 
i*l  i=l  4 

such  that  &(r)  >  1.  (3.26a) 


If  hr»q+l  then  <y(r)  _<  hf  and  from  Lemma  6  we  once  again  get: 


1  B'  (  i,<?(  i)  )  <  2  C'(i,a(i))  for  every  <3  *  it  , , 

i=l  i=l  q  1 

such  that  <y(r)  _>  1.  (3.26b) 

One  readily  sees  that  (3.26a)  and  (3.26b)  are  equivalent  to 
r  r 

2  B '  (  i ,  <y(  i )  +1 )  <  2  C'(i,a(i)+1)  for  every  <y«iT  .  (3.27) 

i=»l  i=l  q 

Following  the  scheduling  from  c0  to  c1#  we  are  left 
with  the  problem  of  scheduling  the  t'^s  from  c^  to  cq+^. 
The  number  of  distinct  due  times  is  now  q  (note  that  t'.  *  0 
for  every  i  such  that  di  =  c1 ) .  Relabel  the  start  time  cx 
as  c'0  and  the  due  times  c2'*-*'Gq+i  as  c'1,...,c,q.  Define 
b"(i,j),  B"(i,j),  and  C"(i,j)  to  be  the  values  obtained  for 

b,  B,  and  C  when  t'^  is  used  in  place  of  t^  and  c' j  is  used 
in  place  of  c^.  We  immediately  see  that  B"(i,j)  = 

B’(i,j+1),  and  C"(i,j)  =  C'(i,j+1),  for  every  i,  j,  l£i£r, 
0^j<c[.  Substituting  into  (3.27)  yields: 

r  r 

2  B"  (i,<y(i) )  2  C"  (i,<y(  i)  )  for  every  G  *  w  . 

i=l  i=l  q 

It  now  follows  from  the  induction  hypothesis  that  the  t'^s 
can  be  scheduled.  [] 

From  Theorem  1,  it  is  clear  that  by  repeatedly  using 
COMPUTE_W  to  determine  the  amount  to  be  scheduled  in  each 
interval,  a  feasible  schedule  can  be  obtained  whenever  such 
a  schedule  exists.  Each  time  COMPUTE_W  is  used,  we  need  to 
recompute  b,  a,  Z,  X,  and  Y.  The  time  needed  for  this  is 
O(nk)(note  that  recurrences  3.6  and  3.10  will  be  used  to 
compute  Z  and  X).  The  for  loop  may  be  executed  in  0(kr  +  n) 
time.  We  may  assume  that  r  _<  n  and  so  the  complexity  of 
COMPUTE  W  is  0(kn) .  The  Kafura-Shen  algorithm  is  of 
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complexity  0(n) .  Hence,  the  overall  computing  time  for  the 
k  phases  of  our  scheduling  algorithm  is  0(k  n) .  An  addi¬ 
tional  O(nlogn)  time  is  needed  to  sort  the  jobs  by  memory 

size  m. .  Hence,  the  overall  complexity  of  our  preemptive 

1  2 
scheduling  algorithm  is  0(k  n  +  nlogn) . 


4. 


Minimizing  Lm 


Let  S  be  a  preemptive  schedule  for  ( t^,d^,m^) ,  1  <_  i  <_  n. 
Let  f^  be  the  finish  time  of  job  i  in  S.  If 

f^  <_  d^,  1  <  i  <  n  then  S  is  a  feasible  schedule  and  no  job 

is  late.  The  lateness  of  job  i  is  f^-d^  and 

Lmax  =  :  l<i<n}  .  Note  that  Ljnax  <_  0  iff  all  jobs 

finish  by  their  due  times.  Also,  note  that  if  L  <  0  then 

c0"cl  i  Lmax ‘ 


From  the  definition  of  Lmax,  if  follows  that  by  chang¬ 
ing  the  release  time  from  c0  to  c0-Lmax  we  obtain  a  job  set 

that  can  be  scheduled  such  that  no  job  finishes  after  its 

due  time.  Hence,  to  determine  the  minimum  L  ,  we  need  to 

determine  the  least  x  such  that  the  condition  of  Theorem  1 

is  satisfied  when  a  release  time  of  c0-x  is  used.  This  x 

may  be  obtained  from  a  form  equivalent  to  that  of  Theorem  1. 

r  r 

We  observe  that  2  B(i,<y(i))^  2  C(i,<y(i))  for  every  <3  *  w  ' 

i=l  i=l 

r  r 

iff  max{  I  B(i,o,(i))-  2  C(i,a(i))}  £  0.  It  is  helpful  to 
o*ir  i=l  i=l 

rewrite  this  form  seperating  out  the  case  when  <3(  i)=0,  1  £  i 

r  r 

<  r.  For  this  <y,  we  see  that  2  B(i,a(i))  -  2  C(i,cr(i))  = 

i=l  i=l 

n 

2b(i,0).  For  every  other  a,  there  is  an  s,  1  _<  s  <_  r  such 
1 

that  <3  (s)  >  1. 


Define  H  as  below: 
s 

s  s 

H  *  max  {  2  B(i,<y(i))-  2  C(i,cr(i))},  l£s£r. 
8  o«n,cr(s)>l  i»l  i*l 
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We  immediately  see  that 


r  r 

maxi  2  B(i,cr(i))  -  2  C(i,<y(i')} 

ct^it  i8*!  i*l 


n 

*  maxlEb(i,0),  H1#  Hr) 


Let  x,  =  max  It-  -  +  ca)  and  let  x?  **  max 

1liln  l_<i<r 

Ih./m.}. 


Clearly,  if  we  change  the  release  time  to  c0-maxl x^ ,  x2 } 
then  max{  Ib'(i,0),  H ‘ ^ ,  H'2,  •••,  H‘r)  =  0  (the  b' ,  H ’ ^ 

values  are  computed  with  respect  to  the  new  release  time 
c0-max{  x1#  x2 } )  .  Hence,  n\axa4nll  B  '  (  i , <3"(  i )  )  -  E  C  '  (i,<J(i) ) } 
0  and  E  B '  ( i , &( i ) )  £  E  C '  ( i , <y(  i) )  for  every  <3  *  it  .  Moreo¬ 
ver,  x  =  maxlx^,  x2 }  is  the  least  value  of  x  for  which  this 
happens.  Hence, 


(Lmax)min  = 


The  H  s  may  be  computed  in  0(kn)  time  as  follows. 

Define  H1  as  below: 
s 


.  o  » 

H*  »  max  {  E  B(j,<y(j))  -  2  C{j,<y(j))J,  l<s^r,  l^i<k. 

8  (S^IT  j=l  j=l 

or(  s) 

Hence,  H  »  H^,  1  <  s  ^  r.  We  immediately  obtain  the  fol- 

S3, 

lowing  recurrence  for  H^: 


max{B(l,j)  -  C ( 1 , j ) }  if  s  *  1 

j>i 

S  3 

E  B(j,k)  -  2  C(j,k)  if  i  »  k 

j»l  j=*l 

max{Hg_^+B(  s,  i)-C(  s,  i) ,  H*+1  ) 


otherwise . 
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Using  this  recurrence,  all  the  s  may  be  obtained  in 
0(rk)  time  (excluding  the  time  needed  to  determine  the 
b(i,j)s,  B(i,j)s  etc.)*  The  additional  time  needed  to  com¬ 
pute  the  B(i,j)s  and  C(i,j)s  is  0(kn  +  nlogn)  (assuming  n  >_ 
m) .  Hence,  the  minimum  L  may  be  determined  in  0(kn  + 

ITlclX 

nlogn)  time.  Having  determined  the  minimum  L  ,  a  schedule 

max 

having  this  Lmax  value  can  be  obtained  by  chaning  Cg  to 

C0~^max^min  an<^  usin9  the  algorithm  of  Section  3. 

5 .  Conclusions 

We  have  developed  an  0(k  n  +  nlogn)  algorithm  to  obtain  a 

preemptive  schedule  for  n  jobs  (t^,d^,m^),  1  _<  i  £  n  on  m 

processors  with  given  memory  sizes.  This  schedule 

minimizes  L  ^  .  The  minimum  value  of  L__„  can  itself  be 
max  max 

obtained  in  only  0(kn  +  nlogn)  time. 
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